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There is a formal correspondence between the isotropic 3- wave kinetic equation and the rate equa- 
tions for a non-hnear fragmentation-aggregation process. We exploit this correspondence to study 
analytically the time evolution of the wave frequency power spectrum. Specifically, we analyzed a 
3-wave turbulence in which the wave interaction kernel is a constant. We consider both forced and 
decaying turbulence. In the forced case, the scaling function diverges as x"^''^ as expected from 
Kolmogorov-Zakharov theory. In the decaying case, the scaling function exhibits non-trivial, and 
hitherto unexpected, divergence with both algebraic and logarithmic spectral exponents which we 
calculate. This divergence leads to non-trivial decay laws for the total wave action and the number 
of primary waves. All theoretical predictions are verified with high quality numerical simulations of 
the 3-wave kinetic equation. 

PACS numbers: 47.35.-1, 82.20.-w, 94.05. Lk 



Wave turbulence is a theory of the statistical evolu- 
tion of ensembles of weakly nonlinear dispersive waves. 
It has been applied to capillary waves on fluid interfaces, 
gravity waves on the ocean, acoustic turbulence and var- 
ious special limits of plasma and geophysical turbulence. 
(For a review of the theory see for a summary of 
applications see 0.) The key feature is the fact that 
weak nonlinearity permits the consistent derivation [3] of 
a wave kinetic equation describing the time evolution of 
the frequency power spectrum, N^^{t). When sources and 
sinks of energy, widely separated in frequency, are added 
to the wave kinetic equation, it can be shown to have 
exact stationary solutions corresponding to a cascade of 
energy through frequency space from the source to the 
sink. The cascade solution is known as the Kolmogorov- 
Zakharov (K-Z) spectrum; it describes an intrinsically 
non-equilibrium state of the wave field. Everything is 
known about the stationary K-Z spectra, their scaling 
exponents, locality and stability. By contrast, very lit- 
tle is known about the time-dependent solutions of the 
wave kinetic equation. A basic scaling theory of the de- 
velopment of the stationary state in the case of forced 
wave turbulence was provided in although numerical 
investigations [Bl, lij have suggested that there are unex- 
plained dynamical scaling anomalies in many cases. Al- 
most nothing is known about time-dependent solutions in 
the case of decaying turbulence where an initial specrum 
is allowed to decay in the absence of external forcing. In 
this Letter we take the first steps to remedy this. 

The subject of aggregation-fragmentation kinetics, 
having its origins in theoretical chemistry has, at first 
sight, rather little to do with waves or turbulence. This 
field concerns itself with the statistical mechanics of en- 
sembles of particles which aggregate or fragment upon 
contact. The principal quantity of interest is the particle 



size distribution, ni{t), denoting the density of clusters 
of mass i at time t. It satisfies a kinetic equation, which, 
in the case of pure aggregation, is the well-known Smolu- 
chowski coagulation equation 7] . For a review of pure ag- 
gregation processes see Q . If clusters also break up, addi- 
tional terms may be added to the Smoluchowski equation 
to take this into account. See [§] for a review of fragmen- 
tation. In aggregation-fragmentation kinetics, in strong 
contrast with wave kinetics, almost all theoretical effort 
has historically been focused on determining the time 
evolution of ni{t) from the underlying kinetic equation. 
As a result, a comprehensive scaling theory of the solu- 
tions of the Smoluchowski equation has been constructed 
(see [13] for a review). Although there is a conceptual 
analogy 11| between energy transfer between scales in 
turbulence and mass transfer between clusters in aggre- 
gation, it is only recently that this analogy has been made 
quantitatively useful. Concepts and techniques from tur- 
bulence have proven useful in analysing certain aspects 



of aggregation problems 
has been shown Il5| that 
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13, Q. 



Furthermore, it 



in the case of isotropic wave 
turbulence with quadratic nonlinearity, the wave kinetic 
equation can be rewritten as a set of rate equations for a 
aggregation-fragmentation process with an unusual non- 
linear fragmentation mechanism. This correspondence 
opens the door for the transfer of ideas and techniques 
from aggregation-fragmentation kinetics to the context 
of wave turbulence which will hopefully start to fill in the 
gap in knowledge of time-dependent solutions of the wave 
kinetic equations alluded to already. Furthermore, this 
correspondence opens up a new set of problems within 
aggregation-fragmentation kinetics. This Letter contains 
some opening explorations in this direction. 

It was shown in [l5| that resonant interactions between 
waves lead to forward transfer of energy between frequen- 
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cies which looks like an aggregation process: (i) © (j) 
{i + j). Back-scatter of energy leads to a fragmen- 
tation process («) (i + j) [i) ® (i) ® (j). This 
fragmentation mechanism is unusual. It is non-linear 
while typically [9] the fragmentation mechanism is linear: 
{i + j) {i) [j). Nonlinear collision-controlled frag- 
mentation processes have been studied in the past (see 
[lit ^^'^ references therein). While they are somewhat 
similar to the above rule, this model asserts that only 
the larger particle breaks and it happens according to a 
rather special rule. Our goal here is to apply ideas and 
techniques developed in studies of aggregation and frag- 
mentation to wave turbulence. We want to examine fun- 
damental aspects and have, therefore, limited ourselves 
to the simplest possible situation where the wave inter- 
action kernels are constant. The dynamical problem is 
already non-trivial at this level. The resulting kinetic 
equation can be reduced (see [15| for details), in the dis- 
crete case llSl, to: 



dt 



W1>1 



where J is the energy injection rate and is the fre- 
quency space wave action. The total wave action is 
N{t) ^ X]td>i ^w(^)- III the decay case (J = 0), it satis- 
fies the equation (found by summing Eqs. p]l) 



dN 
~dt 



(2) 



The exponent r is the wave spectrum exponent or poly- 
dispersity exponent. The shape of the frequency power 
spectrum for large time is determined by the small x 
behaviour of the scaling function, F{x). In aggregation 
problems this divergence has been often encountered and 
the T has proven to be difficult to determine 1^, 17 1: in 
some seemingly simple models the exponent r remains 
unknown. An important lesson from this work is that 
one should be particularly careful when r > 1. 

This is a finite capacity system so there is no dissi- 
pative anomaly / gelation transition. Energy is there- 
fore conserved for all time by the wave interactions. For 
the forced case, the total energy then grows linearly in 
time since we are injecting energy at a constant rate. 
Thus Ml ^ t (we take J = 1). The scaling hypothesis, 
Eq. (HI) , then implies that a = — | and subsequently solv- 

for some constant. 



ing Eq. (O suggests that s{t) ^ 
vq. We then expect the scaling 
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N^{t) ^ s-^/"^ F{uj/s) withs- 



(8) 



On the other hand, for the decaying turbulence energy 
is conserved , Mi(t) — 1. The scaling hypothesis then 
implies that a = —2 and solving Eq. ([5]) gives s{t) ^ sot 
for some constant, sq. We then expect the scaling 



N^{t) - s^^ F{lu/s) with 



Sot. 



(9) 



These predictions for the growth rate of s{t), which are 
based solely on the assumption of scaling and the ab- 
sence of a dissipative anomaly, are verfied numerically in 
Fig. [TJ All numerics have been done using the algorithm 
described in |J5|]. 



The primary waves (monomers) evolve according to 
dNi 

ij>i 
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dt 



uj+l- 



(3) 



We assume the scaling hypothesis: there exists a typical 
scale, s{t), such that N^{t) is asymtotically of the form 

N^it) = s^'Fito/s). (4) 

Given this hypothesis, it follows from Eq. ([1]) that 

ds 



dt 



(5) 



while F{x) must satisfy a complicated integro-differential 
equation. The scale s{t) is defined as a ratio of moments: 



w''N^{t)duj (6) 



Miit) JO 
Often the scaling function, F{x), diverges at small x: 

F{x)^Ax-'' aax^O. (7) 
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FIG. 1: Time evolution of typical frequency, Eq. (|6]), for forced 
and decaying turbulence. The dashed lines correspond to the 
theoretical predictions provided by Eq. ([8]) and Eq. Q re- 
spectively 

For forced turbulence, we expect that the frequency 
spectrum should become stationary for large times. This 
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is an additional piece of information which allows us to fix 
the spectral exponent. Requiring that Eq. ([8]) is indepen- 
dent of t for small u; selects r = 3/2. This corresponds 
to the K-Z exponent for this model [l^. Furthermore, 
the corresponding K-Z constant can be computed exactly 
for this model [15] so that we we obtain asymptotic be- 
haviour of the scaling function 
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2 Vtt - 4 In 2 



as X 



0. 



(10) 



This prediction, and the scaling behaviour, Eq. ([8]) are 
verified explicitly from the numerical data in Fig. [21 




FIG. 2: Time evolution of wave spectrum in the forced case. 
The main panel shows snapshots of Nui{t) at a succession of 
times. The inset shows the same data collapsed according to 
the scaling in Eq. ([8]). The collapsed data has been compen- 
sated by accordance with Eq. (|10p and shows a plateau 
with the theoretically predicted amplitude. 

Let us now turn to the decay case. Figure [3] presents 
numerical simulations of the decay of a monochromatic 
initial spectrum (0) = S^i and verifies the scaling be- 
haviour expected from Eq. The essential difference 
from the forced case is that we no longer have the addi- 
tional constraint provided by stationarity which allowed 
us to easily determine the spectral exponent r. We must 
return to the original kinetic equation. 

Guided by our result for the forced cased, let us pre- 
sume that the wave spectrum diverges as a; — *■ in accor- 
dance with Eq. dH). Substitution of Eq. Q into Eq, ^ 
and comparing requires us to choose: r — 3 = 2t — 4 and 
A = -—^ — 1 This seems to straightforwardly 

determine the spectral exponent to be r = 1 until we re- 
alise that X]w>i [^(J+i) ] ~ 1 resulting in the divergence 
of the amplitude. A, for this choice of r. This surprising 
result suggests that we consider the more general diver- 
gence 



F{x) ~ [ln(l/a;)]^ 



0. 



(11) 
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FIG. 3: Time evolution of wave spectrum in the decay case. 
The main panel shows snapshots of A'^(t) at a succession of 
times. The inset shows same data collapsed according to the 
scaling in Eq. 



where we have introduced a logarithmic spectral expo- 
nent, p, with which one may hope to cancel the diver- 
gence we have just encountered. The tail of the wave 
spectrum then has the form 



A 1 



In 



sit) 



fortj<s(t). (12) 



Setting fc = 1 in this formula gives the asymptotic form 
of of ni{t). Substituting these formulae into Eq. ^ 
one finds that the leading term on the left hand side 
is of order s(t)^^ In [s(t)]'' and the leading order term 
on the right hand side is of order s(i)~^ In [s(i)]^''~^ inot 
sit)~^ In [s(i)]^'' as one might naively expect owing to the 
cancellation alluded to above). Thus we should choose 
p = 1 for the logarithmic spectral exponent so that the 
asymptotic form of the scaling function in the decay case 
is: 
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ln(l/a;) as a; ^ 0. 



(13) 



Fig. m shows the numerically obtained scaling function 
rescaled according to this formula. The plateau at small 
X provides strong numerical support for Eq. p3p . 

In principle, one should also obtain the amplitude at 
this point but this turns out to be easier using Eq. ([2|) 
for Nit) since certain sums which arise can by computed 
exactly in that case. The total wave action is 
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FIG. 4: Scaling function, F{x), compensated for the theoret- 
ical small X divergence, ln(l/a;), expected from Eq. ((TSj. 



Substituting this into the left-hand side of Eq. 0, and 
Eq. (|12p into the right-hand side, and computing the 
leading terms we find the balance 



A so In [s{t)Y 
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The sum gives 7r^/6 from which we conclude that ^ = 

Recalling that s{t) ^ sgt this relation gives us the 
following nontrivial asymptotic decay laws for the total 
density and number of primary waves respectively: 
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(14) 
(15) 



These predictions are validated numerically in Fig. [5] 

To conclude, we have used the analogy between three- 
wave turbulence and aggregation-fragmentation pro- 
cesses to study analytically the decay kinetics of a simple 
wave turbulence model. We found that the kinetics have 
non-trivial scaling properties, even in this simple case, 
which differ significantly from the corresponding aggre- 
gation process. Our results suggest that decaying wave 
turbulence should be studied in greater detail than it has 
been to date. 
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FIG. 5: Time decay of total wave action, N{t), and pri- 
mary wave action, Ni{t), compensated by the theoretically 
predicted asymptotic decay rates taken from Eq. (|14|l and 
Eg. psp respectively. The theoretically expected plateau val- 
ues are attained to within a few percent accuracy. 
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